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■ Abstract 



Influence of weak uniaxial small-scale anisotropy on the stability of inertial-range scaling regimes 
in a model of a passive transverse vector field advected by an incompressible turbulent flow is in- 
U: vestigated by means of the field theoretic renormalization group. Weak anisotropy means that 

pH i parameters which describe anisotropy are chosen to be close to zero, therefore in all expressions it 

is enough to leave only linear terms in anisotropy parameters. Turbulent fluctuations of the velocity 
field are taken to have the Gaussian statistics with zero mean and defined noise with finite correla- 
tions in time. It is shown that stability of the inertial-range scaling regimes in the three-dimensional 
case is not destroyed by anisotropy but the corresponding stability of the two-dimensional system 
^ I can be destroyed even by the presence of weak anisotropy. A borderline dimension dc below which 

■ the stability of the scaling regime is not present is calculated as a function of anisotropy parameters. 

O . 

^ Introduction 

O I During the last decade models of scalar or vector fields passively advected by prescribed stochastic 

■ environment have played the central role in the process of understanding of the so-called anomalous 
scaling, the term that refers to the possible deviations from the predictions of the Kolmogorov phe- 
nomenological theory [H [2] . Theoretical understanding of the anomalous scaling in the framework of a 

^ ^ ■ microscopic model remains one of the main unsolved problems in the theory of fully developed turbu- 
^ . lence. On the other hand, it is well known that the breakdown of the classical Kolmogorov-Obuchov 
^ I phenomenological theory of fully developed turbulence [2] is even more noticeable for simpler models 
of passively advected scalar or vector quantity than for the velocity field itself and, at the same time, 
the problem of a passive advection is easier from theoretical point of view (see, e.g., [3] and references 
therein). 

One of the most suitable approach for studying self-similar scaling behavior is the method of 
the field theoretic renormalization group (RG) [H [5] which was widely used in the theory of critical 
phenomena. This method can be also used in the theory of fully developed turbulence and related 
problems [5l[6l[7], e.g., in the problem of a passive scalar (or vector) field advected by a given stochastic 
environment. In [8] the field theoretic RG was first time applied to the model of a passive scalar 
advected by a given statistics of velocity field, namely, to the so-called Kraichnan model [9j, where 
a scalar field is advected by a self-similar white-in-time velocity field. It was shown that within the 
field theoretic RG approach the anomalous scaling is related to the existence of "dangerous" composite 
operators with negative critical dimensions in the framework of the operator product expansion (OPE) 
[5l[6l[7]. Afterwards, various generalized descendants of the Kraichnan model, namely, models with 
inclusion of large and small scale anisotropy, compressibility, and finite correlation time of the velocity 
field were studied by the field theoretic approach (see [TO] and references therein). Moreover, advection 
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of a passive vector field by the Gaussian self-similar velocity field (with and without large and small 
scale anisotropy, pressure, compressibility, and finite correlation time) has been also investigated and 
all possible asymptotic scaling regimes and crossovers among them have been classified |lllll2^ [T3 l ll4j. 
General conclusion is: the anomalous scaling, which is the most important feature of the Kraichnan 
rapid-change model, remains valid for all generalized models. 

In [15] the influence of the small-scale anisotropy on the infrared (IR) stability of the possible 
scaling regimes was investigated in one particular model of a passive vector advected by a Gaussian 
velocity field with finite correlation time in the presence of the small-scale anisotropy, namely the 
model where the stretching term is absent (the so-called ^ = model, see, e.g, [121 114]). It can be 
consider as a starting point for studying of the influence of anisotropy on the anomalous scaling of the 
model. But, as for anomalous scaling, this model is a little bit specific because in contrast to the other 
models of passive vector admixture where the anomalous scaling is related to the composite operators 
built of the vector field without derivatives |13[ [T3] in the case under consideration it is related to 
the composite operators built solely of the gradients of the field. This fact radically changes the 
complexity of the problem especially in the anisotropic case (see, e.g., \12\ 116] and references therein). 
Thus, as it was stressed in |15] . it can be consider as a further step to the nonlinear Navier-Stokes 
equation problem. Because the problem is complicated even at the first stage of the RG analysis (see 
[15]). in what follows, we shall return to the problem of the influence of the small-scale anisotropy on 
the IR scaling regimes, namely, we shall try to understand it when the anisotropy is considered to be 
weak. It means that all nonlinear terms in respect to anisotropy parameters in all expressions can 
be neglected. In this case, one has explicit expressions for all quantities (coordinates of fixed points, 
eigenvalues of a matrix of the first derivatives, etc.) and the analysis of the possible scaling regimes 
can be done analytically. The results of the present paper will be used in the subsequent investigations 
of the anomalous scaling in the model. 



Formulation of the model 

The model of the advection of transverse (solenoidal) passive vector field b = b(x, t) is described by 
the following stochastic equation 

dth = z^oAb - (v • V)b + f , (1) 

where dt = d/dt, A = is the Laplace operator, vq is the diffusivity (a subscript denotes bare 
parameters of unrenormalized theory), and v = v(x, t) is incompressible advecting velocity field. The 
vector field f = f (x, t) is a transverse Gaussian random (stirring) force with zero mean and covariance 

dI = (/,(x,t)/,(x',t')) = Sit - t')a,ir/L), r = X - x' (2) 

where parentheses (...) hereafter denote average over corresponding statistical ensemble. The noise 
defined in ([2]) maintains the steady-state of the system but the concrete form of the correlator will not 
be essential in what follows. The only condition which must be satisfied by the function Cij{r/L) is 
that it must decrease rapidly for r = |r| ^ L, where L denotes an integral scale related to the stirring. 

In real problems the velocity field v(x) satisfies Navier-Stokes equation but, in what follows, we 
shall work with a simplified model where we suppose that the velocity field obeys a Gaussian statistics 
with zero mean and pair correlation function 

{v.{x)v,{x')) ^ D^^{x;x') = J i^i2,,.(k)i?"(^,k)e--(*-*')+^'^(---'), (3) 

where d is the dimension of the space, k is the wave vector, and Rij{\s.) is a transverse projector. In 
our uniaxial anisotropic case it is taken as (see, e.g., [13] and references therein) 

Rijik) = (l + Qi(n • kf/k^) P,j(k) + a2nsniPis(k)Pji{k) , (4) 
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where -Pj/;(k) = Sij — kikj/k"^ is common isotropic transverse projector, the unit vector n determines the 
distinguished direction, and qi, 02 are parameters characterizing the anisotropy. From the positiveness 
of the correlation tensor D^j one finds restrictions on the values of the above parameters, namely, 
ai^2 > —1- The function Z'''(a;,k) in ^ is taken in the following form [14j 

D^iu k) = 9oUouik^-'~'^-v 

^ ' ' {iuj + UQVok'^-^)(-iuj + uou^k'^-'^y ^ ' 

where qq plays the role of the coupling constant of the model (a formal small parameter of the ordinary 
perturbation theory), the parameter uq gives the ratio of turnover time of scalar field and velocity 
correlation time, and the positive exponents e and r] are small RG expansion parameters. The coupling 
constant and the exponent e control the behavior of the equal-time pair correlation function of 
velocity field and the parameter uq together with the second exponent r] are related to the frequency 
LO ~ UQVQk'^~^ which characterizes the mode k. The value e = 4/3 corresponds to the celebrated 
Kolmogorov "two-thirds law" for the spatial statistics of velocity field, and r/ = 4/3 corresponds to 
the Kolmogorov frequency. Dimensional analysis shows that and uq, which we commonly term as 
charges, are related to the characteristic ultraviolet (UV) momentum scale A (or inner length / ~ A^^) 
by 

50 ^ A2^ uo ~ A^. (6) 
The stochastic problem ([IJ-Q can be rewritten in a field theoretic form with action functional 

mm 

S{^) = \)j (-(9t - WjcJi i^oA + i^oXio(n • 9)^) + "-i t'o (x2oA + X3o(n • 9)^) nfc 6^. 

-\{v,\DlY\,-U,Dy^. (7) 

where and are given in ^ and ([2]) respectively, b' is an auxiliary vector field (see, e.g., [5]), 
and the required integrations over x = (x, t) and summations over the vector indices are implied. In 
action ([7]) the terms with new parameters Xio 1X205 and xso are related to the presence of small-scale 
anisotropy and they are necessary to make the model multiplicatively renormalizable [5]. Model ([7]) 
corresponds to a standard Feynman diagrammatic technique (see, e.g., |13j for details) and the stan- 
dard analysis of canonical dimensions then shows which one-irreducible Green functions can possess 
UV superficial divergences. Detail analysis of the RG technique in the model will be given elsewhere. 
We stress only that the functional formulation ([7]) gives possibility to extract large-scale asymptotic 
behavior of the correlation functions after an appropriate renormalization procedure which is needed 
to remove UV-divergences. 



Influence of anisotropy on scaling regimes of the model 

The RG analysis leads to the conclusion that possible scaling regimes are given by the IR stable fixed 
points of the corresponding RG equations O [6l [?]. The fixed points of the RG equations can be 
determined in two ways: First, they are determined by the corresponding system of RG differential 
equations (they are know as the flow equations or Gell-Mann-Low equations), or, second, they can 
be found from the requirement that all the so-called beta functions of the model vanish and the IR 
stability of the fixed point is given by the requirement that all the eigenvalues of the matrix of the 
first derivatives Vlij = d(3i/dCk must have positive real parts, where (3i denotes the full set of beta 
functions and is the full set of charges of the model. In the case when uniaxial anisotropy is 
strong (nonrestricted), only the first possibility is suitable. It was briefly discussed in [15j where the 
corresponding analysis of the scaling regimes was done. In what follows, we shall concentrate on the 
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Figure 1: (Left) The scaling regimes of the model in the e — rj plane. The region FPI corresponds to 
the trivial rapid-change limit: g* /u* = 0, 1/u* = 0. The region FPU corresponds to the the non-trivial 
rapid change limit: g* /u* > 0, 1/u* = 0. The region FPIII corresponds the trivial "frozen" limit: g* = 
0,u* = 0. The region FPIV is related to the non-trivial "frozen" limit: g* > 0,u* = 0. In the end, the 
line £ = r] (FPV) corresponds to the more interesting scaling regime with g* > 0,0 < u* < oo. (Right) 
Dependence of the borderline dimension dc on weak anisotropy parameters ai^2 in the Kraichnan limit 
of the model: oo. 



weak anisotropy limit to better understand the situation. In this case, the second possibility leads to 
the result, i.e., the possible fixed points are given by the following system of equations 

Pg{g*,u*,Xi*,aj,d,£,r]) = [3u{g*,u^,Xi*,(y-j,d,e,vi) = (3^^{g^,u^,Xi*,aj,d,£,ri) = 0, (8) 

for i = 1, 2, 3; J = 1, 2, where each variable with the star denotes the corresponding fixed point value of 
the variable and all beta functions are linear functions of all quantities related to anisotropy, namely, 
ai,a2, and Xii ^ = li 2, 3. The explicit form of the beta functions is as follows 

f3g = gi-2e + 2ji), /3„ = «(-?? + 71), /^^^ = Xi(7i - 7m), ^ = 1,2,3, (9) 

where, in the limit of weak anisotropy, the so-called gamma functions 7^, i = 1,2,3,4 have the following 
explicit form 

= "2^ {27Tr d{d + 2){d + 4)(d + 6)(1 + ' ^ = 0'1'2'3 (10) 

where we define xo = = 2'k'^/'^ /T{d/2) is the d— dimensional sphere, and the coefficients /Ci,i = 
1, 2, 3, 4 are given as follows 

JCi = -(02(3 + (i)(6 + d){l + u) + ai(6 + d){l + d{A + d)){l + u)- 3(24 + 2xi + 6x2 + 5^3 + 24n) 
+d(-30 - 25x1 - 9x2 - 3X3 - 30n + d(21 + lOd + d^ - (10 + d)xi - X2 + (3 + d)(7 + d)u))) 

JC2 = -2(6 + (i)((3 + d)xi + X2) + 12x3-2ai(3 + (i)(6 + d)(l + u) 
+02(6 + d){-lO + d{-6 + d{3 + d)))(l + u) , 

IC3 = (6 + d){-2{3 + d)xi + i-W + d{-6 + d{3 + d)))x2) 
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Figure 2: Dependence of the borderline dimension dc on weak anisotropy parameters ai_2 in the case 
of the so-cahed frozen hmit of the modeh = 0, and in the case with finite value: = 1. 

+(-24 + d{-8 + d{5 + d)))x3 + 2a2(6 + d){l + u) + 2ai(3 + d){6 + d){l + u) , 
/C4 = -2d((i + 4)x3. 

Thus, possible fixed points are found as solutions of the system of algebraic equations ([8]) and their 
IR stability is determined by the positive real parts of the eigenvalues of the matrix Q = {Qik}- First 
of all, we have found all possible IR fixed points which correspond to the possible scaling regimes and 
we have analyzed the regions of their IR stability in the e — t] plane. The results of this analysis are 
shown in Fig.[T]||Left), where it is shown that the model exhibits five diff'erent scaling regimes (two for 
rapid-change limit, two for so-called "frozen" limit, and one general with nonzero u*) (see, e.g., jl4J and 
references therein). The same situation is also held when no restrictions on the uniaxial anisotropy is 
supposed [15]. Further, our interest is concentrated on the investigation of the dependence of stability 
of the above mentioned scaling regimes on anisotropy parameters 01,02 and on the dimension of the 
space d. As in the strong anisotropy case we have found the so-called borderline dimension dc between 
stable and unstable regimes as a function of anisotropy parameters 01,02 and parameter u^,. The 
results are shown in Figs.l^Right) and [2] for different fixed point values of the parameter u. The 
assumption of linearity of beta functions as a functions of anisotropy parameters leads to the fact 
that the borderline surface dc = dc{cti,Oi2) consists of two intersecting planes. One of them is related 
to the condition > and the second surface is related to the condition to have all real parts of 
eigenvalues of the matrix of the first derivatives positive. One can see that the presence of small-scale 
anisotropy leads to the violation of the stability of the corresponding scaling regimes below dc € [2, 3] 
for appropriate values of anisotropy parameters. But from the point of view of further investigation 
of anomalous scaling the most important conclusion is that all the three-dimensional scaling regimes 
remain stable under influence of small-scale uniaxial anisotropy. 

Conclusions 

In present paper we have studied the influence of weak uniaxial small-scale anisotropy on the stability 
of the scaling regimes in the model of a passive vector advected by given stochastic environment with 
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finite time correlations by means of the field theoretic RG. The system exhibits five different scaling 
regimes. Their IR stability is related to the values of the parameters e and rj. Besides, the stability 
of all scaling regimes is influenced by presence of small-scale anisotropy which is demonstrated in the 
existence of the so-called borderline dimension dc which is a function of the anisotropy parameters. 
The dc is defined as dimension above which the corresponding scaling regime is stable and below which 
the stability of the regime is destroyed. We have calculated the borderline dimension in the case when 
anisotropy parameters are close to zero. The assumption of linearity of beta functions as a functions 
of anisotropy parameters leads to the fact that the borderline surface dc = dc{ai,a2) consists of two 
intersecting planes. All the calculations have been done at the one-loop level. The results will be used 
in the further investigations of the anomalous scaling of the model. 
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